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ABSTRACT. - The species Candidia barbata , collected from different river basins in Taiwan, was surveyed 
for mitochondrial genetic haplotypes. The species is endemic to Taiwan and is widely-distributed around the 
island. A total of 208 samples of C. barbata specimens yielded 53 mitochondrial haplotypes of D-loop 
sequences. Using congruent tree topology, two major clades were detected from Taiwan: Clade A (North and 
West) and Clade B (South). The results were obtained by employing both distance (neighbor-joining) and 
discrete (maximum parsimony) methods of sequence analysis. The mitogenetic diversity is high in this species 
complex. 

KEY WORDS. - Candidia barbata, Cyprinidae, species complex, molecular evolution, biogeography, Taiwan. 


INTRODUCTION 

The small-sized cyprinid genus Candidia (Jordan & 
Richardson, 1909) has been regarded as endemic to Taiwan 
(Chen, 1982; Tseng, 1986; Chen & Fang, 1999). However, 
it shares close resemblance with the Japanese Zacco taxon. 
The species Candidia barbata (Regan, 1908) was originally 
described as Opsariichthys barbatus. However, the species 
was validated in the genus Candidia by Chen (1998). Chen 
& Fang (1999) listed the distribution pattern of Candidia in 
Taiwan and it can be found from Taipei County to Pingtung 
County, mainly in the Western drainages next to the Central 
Mountain Ridge of Taiwan. This species inhabits small 
tributaries of river basins (Chen & Fang, 1999; Chen & 
Chang, 2005). 

Taiwan, the largest island off the mainland of China, was 
formed and lifted by the interaction of the Euro-Asia and 
Philippine Sea Plates. The waters are shallow on the 
Formosan bank on the West and very deep (up to 1,000 - 
2,000 m) on the Eastern slope. The freshwater river systems 
mainly flow to the Western and Eastern coasts, being well 
separated by the Central Mountain Ridge (Page & Suppe, 
1981). 


The Order Cypriniformes, being primarily freshwater fishes, 
are strongly affected by the glaciations and river connection 
via land bridges into the various river basins between 
mainland China and the island of Taiwan. Low temperatures 
affect the inland-water fish communities, causing the 
communities to migrate to the lower reaches from the 
mainland river systems to avoid temperature stress. The 
communities then connect to nearby tributaries of other rivers 
which were once joined as ancient river basins (Froufe et al., 
2002; Kotlik & Berrebi, 2002; Perdices et ah, 2003). 

The mitochondrial DNA (mtDNA) exhibits a high 
evolutionary rate with the accumulation of neutral mutations 
and thus, accounts for evolutionary events. Herein, we use 
the D-loop sequence (or control region) as the DNA genetic 
marker. This sequence has an AT-rich non-coding region in 
most species. Since the sensitivity of such DNA markers is 
good enough to detect the polymorphic haplotypes among 
species, it is possible to reconstruct the distribution pattern 
of wide-ranging species by using more detailed molecular 
phylogeographic and morphological studies (Chen et ah, 
1998, 2002). This paper reports on the genetic diversity and 
geographical speciation patterns of the C. barbata species 
complex in Taiwan and hypothesizes as to whether the 
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speciation event of this species complex is based on 
monophyletic species or resulting from the isolation of sibling 
species. 

MATERIALS AND METHODS 

A total of 208 individuals of Candidia barbata were collected 
from 29 river basins in Taiwan (Fig. 1 & Table 1). Two 
morphologically similar species, Zacco temminckii 
(Temminck & Schlegel, 1846) and Zacco sieboldii 
(Temminck & Schlegel, 1846), which have been recently 
revalidated by Hosoya et al. (2003), were assigned as 
comparative outgroup taxa. 

All DNA samples were extracted using the general protocols 
of the phenol-chloroform method (Sambrook et al., 1989; 
Chen et al., 2002). The DNA fragments of about 1,500 base 
pairs (bp), inclusive of the full length D-loop region were 
amplified by polymerase chain reaction (PCR) using a set of 
primers for the flanking regions which were designed from 
the sequences of tRNA-PHE and 12S rRNA (Chen et al.. 
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Fig. 1. Sampling sites of Candidia barbata for mitogenetic D-loop 
sequences analysis, showing Grades I to VI. The gray regions depict 
mountainous areas that might be potential geographical barriers. 
Numbers represent rivers and names are indicated in Table 1. 


1998, 2002). PCR was done in a Model 9700 thermal cycler 
(Perkin-Elmer) with 30 - 40 cycles carried out for each 25 
uL reaction volume which contained 14.4 iiL sterile distilled 
water, 2.5 uL lOx PCR buffer (Takara), 2.0 uL dNTP (2.5 
pM each), 2.5 uL each primer (5 uM each), 0.1 pL 0.5 units 
Ex Taq (Takara) and 1.0 pL template DNA. Each thermal 
cycle profile was as follows: denaturation at 94°C for 15 
seconds, annealing at 50 °C for 15 seconds and extension at 
72°C for 60 seconds. The PCR products were electrophoresed 
on a 1.0% agarose gel (Takara) and which was then stained 
with ethidium bromide for band characterization via 
ultraviolet trans-illumination. 

The double-stranded PCR products were purified using a kit 
(Boehringer Mannheim, High Pure PCR Product Purification 
kit) and were subsequently used for direct cycle sequencing 
with florescent dye-labeled terminators (ABI Big-Dye kit). 
The primers used were the same as those used for PCR. All 
sequencing reactions were performed according to the 
manufacturers’ instructions. Labeled fragments were analyzed 
using an ABI Model 377-64 DNA sequencer (ABI). 

Distinction of divergent haloptypes of all nucleotide 
sequences for Candidia species were conducted using 
Molecular Evolutionary Genetics Analysis program (MEGA) 
version 3.0 (Kumar et al., 2004). Evolution model testing was 
conducted using Modeltest version 3.7 (Posada & Crandall, 
1998). Phylogenetic and molecular evolutionary analyses 
were performed using PAUP version 4.0b 10 (Swofford, 
2002 ). 

The minimum-spanning network was constructed with the 
MINSPNET (Excoffier & Smouse, 1994). Inter- and intra¬ 
population genetic diversity were quantified by indices of 
haplotype diversity (Nei & Tajima, 1983) and estimates of 
nucleotide diversity (Jukes & Cantor, 1969) using DnaSP 
version 3.14 (Rozas & Rozas, 1999). 


RESULTS 

A total of 208 Candidia barbata samples from 29 different 
river basins in Taiwan were analyzed with two outgroups 
(Zacco temminckii and Z. sieboldii). There were 53 distinct 
haplotypes from the mtDNA D-loop sequences which could 
be divided into two distinct clades. Clade A comprises of 42 
halotypes (GenBank accession numbers: DQ489379 to 
DQ489420) and Clade B comprises of 11 haplotypes 
(GenBank accession numbers: DQ489421 to DQ489431). 

The length of complete D-loop sequences is 995 - 997 bp in 
Clade A, 996 bp in Clade B, 996 bp in Zacco temminckii and 
997 bp in Z sieboldii. The aligned mtDNA sequences for all 
haplotypes were up to 1,000 bp. The most common 
haplotypes were as follows: All, A30 (both based on 32 
individuals), A26 (based on 25 individuals) and A40 (based 
on 16 individuals). The other haplotypes each occurred in 
less than 10 individuals (Table 1). Among the distribution 
patterns of haplotypes, the most widely-distributed was Al 1 
(found in seven river basins), followed by A30 (found in four 
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Table 1. Distribution of haplotypes from mitochondrial D-loop sequence of Candidia barbata complex in Taiwan. 
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river basins), A12 and B06 (each found in three river basins) 
and the rest were found in less than two river basins. 

The phylogenetic analysis, reconstructed from the different 
haplotypes, was conducted using the distance method 
[neighbor-joining (NJ) method based on the GTR model 
which was chosen by Modeltest from the sequence data 
matrix] and the discrete method [maximum parsimony (MP) 
method by gap complete deletion with 989 bp]. Both methods 
yielded very similar topology and supported the same six 
grades (Clade A: Grade I, II, III, IV, V; Clade B: Grade VI) 
(Fig. 2). The MP tree had 136 trees with equal minimal lengths 
of 279. The majority-rule consensus shows the consensus tree 
with NJ tree topology with 1,000 bootstrap replications (Cl 
= 0.62, R1 = 0.91). The minimum spanning network was also 
reconstructed by the six distinct grades and a long branch for 
Clade B (Fig. 3). From the estimation for divergence times 
of mitochondrial D-loop sequences by Johns & Avise (1998), 
this was 5% per million years. Cortey & Garcia-Marin (2002) 
estimated that the Atlantic brown trout divergence rate was 
2%. Wang et al. (2004) used the molecular clock for the 
cyprinid, Varicorhinus barbatulus (= Onychostoma 
barbcibtulus) at around 3.2 - 10.4%. 


Grade Clade 



"\ 

I 


II 





VI 


B 


Fig. 2. Neighbor-joining tree of mtDNA haplotypes based on 
sequence variation of the D-loop region in Candidia barbata 
showing the divisions into six grades (I - VI). Numbers at the nodes 
are the bootstrap percentage values for 1,000 replicates by neighbour¬ 
joining (NJ) and maximum parsimony (MP) methods. Zacco 
temminckii (ZTE) andZ. sieboldii (ZSI) are used as the outgroups. 


This study employed the estimation as an average of genetic 
divergence of 6% per million years. Hence, the separation of 
Clades A and B can be estimated to have occurred 450,000 
years ago and the separation between all the C. barbata 
populations and the assigned outgroup, the Z. temiinckii 
complex, was 750,000 years ago. All of the events happened 
after late Pleistocene. 


DISCUSSIONS 

The Candidia barbata species complex in Taiwan is more 
similar to the Japanese Zacco temminckii and Z sieblodii in 
overall morphology and genetic homology than to other 
members of Zacco (Hosoya et al., 2003). These groups share 
a similar preference of habitat, inhabiting the upper reaches 
of rivers and small tributaries rather than the main river basins. 
However, congeneric or sibling species of Candidia have yet 
to be found in mainland China. From our unpublished data 
based on recent surveys, Z. temminckii may be considered 
the sister group of the Taiwanese populations of Candidia. 
The generic position of Z. temminckii is still in need of 
clarification. More detailed morphological and osteological 
studies will be needed to decide if Z. temminckii is to be placed 
in a genus separate from Zacco or reassigned as a Candidia 
species. 

The six-grade division was generally congruent with the 
geographic pattern of Taiwan. Clade B (Grade VI) in 
Henchuen Peninsula is well-isolated from the other 



Fig. 3. The minimum-spanning network made by mtDNA 
haplotypes of the Candidia barbata complex. I - VI = grades; 
haplotype (small circle) =1-10 individuals; haplotype (medium 
circle) =11-20 individuals; haplotype (large circle) = more than 
20 individuals. 
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geographical regions. Grade V is limited to Northern Taiwan, 
which is well-separated from the other grades (I, II, III and 
IV) in Clade A. Among them, Grades I and II are much closer 
to each other. Grade IV is within the Southern region of the 
Alisan Ridge which is far from Grade III in the Northern 
region of the Alisan Ridge. This separation of the Southern 
Clade B can be explained by the long-term geographical 
isolation and leading to a greater differentiation between such 
populations, which may be considered as discrete allopatric 
species. 

Analysis of the sediment Oxygen-18 isotope ( ls O) records 
around the South China Sea have revealed the occurrence of 
20 glacial periods (marine isotope stages, MIS) within the 
past 800,000 years (Boulay et al., 2003). From the mitogenetic 
divergence rates of C. barbata, the earliest separation between 
the Taiwanese C. barbata complex and the Japanese Z. 
temminckii and Z. sieboldii can be estimated to be around 
750,000 years (during the MIS 17 period). The separation 
between Clade A and B was around 450,000 years ago (at 
the start of MIS 11), which can be considered the glacial 
period during MIS 12 and allowed connection by land 
between mainland China and Taiwan by sea level regression. 
This land connection might be the isolating mechanism for 
Clades A and B. After this glacial period, the coastline of the 
Southern region of Taiwan seems to be well-isolated between 
these geographical regions during the inter-glacial period and 
all the events of the glacial periods that follow. 

The species complex of two distinct clades was well supported 
by both distance and discrete methods. The mitogenetic 
diversity is still very high suggesting a good genetic pool 
still exists. The different geographical regions seem to have 
independently-evolved mitogenetic haplotypes as Clade B 
was well-differentiated from Clade A from the molecular 
sequence evidence. The allopatric group, Candidia sp. in the 
South, had been recently examined in detail and exhibits 
distinct morphological and meristic features and is currently 
being described (Chen et al., unpublished data). The 
remaining main group, inclusive of the North and West 
groups, would thus be assigned as C. barbata. The occurrence 
of this endemic species complex in Taiwan serves as an 
important example when conservation of freshwater fish 
communities should be considered. 
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